Asymptotic shape of small cells 



Mareen Beermann* Claudia Redenbach^and Christoph Thale* 



Abstract 

A stationary Poisson line tessellation is considered whose directional distribution is concentrated 
on two different atoms with some positive weights. The shape of the typical cell of such a tessellation 
is studied as its area or its perimeter tends to zero. In contrast to known results where the area 
or the perimeter tends to infinity, it is shown that the asymptotic shape of cells having small area 
is degenerate. Again in contrast to the case of large cells, the asymptotic shape of cells with small 
perimeter is indeterminate. The results are accompanied by a large scale simulation study. 
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1 Introduction 

In the preface of the book |15| . D.G. Kendall re-phrased a conjecture about the shape of planar 
tessellation cells having large area. He considered a stationary and isotropic Poisson line tessellation 
in the plane and conjectured that the shape of the cell containing the origin is approximately circular 
if its area is large. First contributions to Kendall's conjecture are due to Goldmann [6j, Kovalenko 
|10| and Miles [12] . The first result for higher dimensions is by Mecke and Osburg [TT] who considered 
so-called Poisson cuboid tessellations. In a series of papers, Calka [I], Calka and Schreiber [3j, Hug, 
Reitzner and Schneider [7j and Hug and Schneider [U [9] treated very general higher-dimensional 
versions and variants of Kendall's problem for quite general tessellation models (Poisson hyperplanes, 
Poisson- Voronoi and Poisson-Delaunay tessellations) and size functionals, see also the book chapter 
[2] for an overview. The respective results either rely on asymptotic theory for high density Boolean 
models or on sharp inequalities of isoperimetric type. 

In this paper, we focus on the analysis of the shape of small tessellation cells. So far, we were 
not able to discover a general principle as mentioned above for the large cells behind the asymptotic 
geometry of small cells. For this reason, we restrict attention to the following simple model and 
its affine images. Take two independent stationary (homogeneous) Poisson point processes of unit 
intensity on the two coordinate axes in the plane and draw vertical lines through the points on the 
x-axis and horizontal lines through the points on the y-axis; see Figure [l] (left). The collection of 
these lines (without the two coordinate axes) decomposes the plane into a countable number of non- 
overlapping rectangles, the collection of which is called a rectangular Poisson line tessellation; see 
[5] . Of interest in the present paper is the shape of a typical rectangle of the tessellation (the precise 
definitions follow below). Mecke and Osburg [11] have shown that a typical rectangle tends to be 'more 
and more cubical as the area tends to infinity'. In the present paper we are interested in the converse 
question and ask for the shape of a typical rectangle of small area. We will show that, in contrast 
to the large area case, the shape of typical rectangles with small area is asymptotically degenerate. 
Besides rectangles of small area, we also consider rectangles that have small perimeter. For such a 
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Figure 1: A rectangular Poisson line tessellation (left) and a non-orthogonal Poisson parallelogram tessellation 
(right). 



situation we obtain a uniform distribution for our parameter measuring the shape of the rectangles, 
leading to an indeterminate asymptotic shape. Again, this result is in contrast to the large perimeter 
case indicated in [TT] and rigorously analyzed in [13]. We would like to stress the fact that this is 
the first paper dealing with the mathematical analysis of small cells in random tessellations although 
some conjectures together with heuristic arguments have appeared earlier in |12] . 

The rest of this text is structured as follows. In Section [2] we present the mathematical framework, 
the statements of our results, Theorem [T] and Theorem [2j a large scale simulation study as well as an 
outlook to higher space dimensions. The proof of Theorem [T] is the content of Section [3] and Theorem 
[2] is derived in Section HI 

2 Results 
2.1 Framework 

Denote by C the space of lines in M? and by Co the subspace containing the lines through the origin. 
We let L\ and Li be two different lines in Co and fix q £ (0, 1). On Co we define the probability 
measure Q by Q = qd^ + (1 — q)^l 2 i where bj Ji stands for the Dirac measure concentrated at Lj, 
% = 1,2. This is to say, Q is concentrated on L\ and Li with weights q and 1 — q, respectively. We 
also define the translation invariant measure A on £ by the relation 



where £ G ± stands for the Lebesgue measure on G- 1 , 7 € (0, 00) and /:£—)• M is a non- negative 
measurable function. In other words, A is concentrated on two families of lines parallel to L\ and L2, 
whereas 7 is an intensity parameter. 

Let now r\ be a Poisson point process on C with intensity measure A as in Q; cf. |14| [T5] for 
definitions. Clearly, the lines of rj decompose the plane into countably many parallelograms - called 
cells in the sequel - which have pairwise no common interior points; see Figure[T] (right). The collection 
of all cells is denoted by C = C{rj) and the intersection point of the two diagonals of a parallelogram 3? 
is denoted by c(^). We may define a probability law P on the (measurable) space of parallelograms 
as follows: 



(1) 




E J2 - c(&>) G ^4} l{c(^) C [n]} 



(2) 



P(^4) = lim 



E £ l{c(^) C [n]} 
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where A is a measurable subset of parallelograms and [n] stands for the centered unit cube of volume 
n. The definition ([2]) formalizes the idea of a uniformly selected cell from C (regardless of size and 
shape) and we call a random parallelogram with distribution P a typical cell of the tessellation; cf. 
|14| [To] for background material. 

2.2 Results for small cells 

We consider the typical cell of a Poisson line tessellation as described above. Its random edge lengths 
are denoted by X and Y and its area by A = XY. To measure the shape of the typical cell we 
introduce two deviation Junctionals. The first one is 

min{X,y| 

(3) " = 2 ^TY 1 ' 

which is a random variable taking values in [0, 1] (this was the reason for the choice for the factor 2). 
We notice that a is scale invariant, i.e., a does not change if the parallelogram is rescaled by some 
constant factor. Moreover, we have a = if at least one of the edge lengths X or Y is zero, i.e., if the 
parallelogram degenerates in that it is a line segment or even a single point. As a second deviation 
functional we introduce 

(4) r = max{X,y}. 

This is not a scale invariant quantity, but we notice that r = if and only if the parallelogram is 
degenerate to a point. 

We investigate first the asymptotic behavior of the deviation functionals a and r under the condi- 
tion that the typical cell area A tends to zero. Our main result in this direction reads as follows. 

Theorem 1. Let < e < 1. It holds that 

(5) P(<7 > e\A <a) = 0{^fa) as a -> 
and consequently 

(6) limP(cj > e\A < a) = 0. 

a— s>0 

Moreover, 

(7) lim P(t > e\A < a) = 0. 

a— >0 

Some comments are in order about the interpretation of Theorem [TJ Firstly, ^ shows that the 
asymptotic shape of a typical cell of small area tends to that of a (possibly degenerate) line segment. 
On the other hand, ([7]) shows that, in the limit, this line segment cannot have positive length. Together, 
Q and ([7]) show that the edge lengths of the typical cell decrease at different rates as A — > 0. In other 
words this means that in a first stage the typical cell shrinks basically in only one direction, which 
causes a typical parallelogram with small (but positive) area to look like a line segment. However, 
in a second stage, as the area shrinks even more, this line segment eventually degenerates to a point. 



This phenomenon is well reflected in the simulation study presented in Section 2.3. We also remark 
that ((5j) gives an upper bound for the rate of convergence. 

Besides cells of small area, also cells with small perimeter can be considered. In this case the 
picture is somewhat different from that presented for the small area case in Theorem [TJ In what 
follows we denote by P half of the perimeter length of the typical cell (the factor 1/2 is chosen for 
simplicity as will become clear in the proof). 

Theorem 2. Let < e < 1. Then a is uniformly distributed on [0, 1] given that P < p, i.e., 

F(a >e\P<p) = 1-e, 
independently of p. Moreover, for any p < e we have that P(r > e\P < p) = so that 
(8) limP(r > e\P <p) = 0. 
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Figure 2: The 150 cells with smallest area (left) and smallest perimeter (right). Note that the axes on the 
left are cut off. The maximal edge length observed is 0.09. 



Theorem [2] shows that the conditional deviation functional a, given P < p, follows a uniform 
distribution on its range [0, 1]. This means that not only in contrast to the case of large perimeter (see 

I13| ) . but also in contrast to the case of small area considered in Theorem [l] above, the asymptotic 
shape of cells that have small perimeter is indeterminate. This phenomenon is well reflected by 



the simulation study in forthcoming Section 2.3 We also refer to a related short discussion at the 
beginning of Section 7 in [9j about the independence of the shape of the typical cell and its perimeter. 
Moreover, Q shows that in contrast to Theorem [I] the deviation functional r does not provide useful 
information about the asymptotic cell shape in the present situation since P(r > e\P < p) = for any 
p < e. 



2.3 Simulation results and outlook to higher space dimensions 

To highlight and to underpin the theoretical results in the previous subsection we performed the 
following simulation study. We choose 7 = 2, q = 1/2 and Q = \(&l x + $L V ), where L x and L y 
are the two orthogonal coordinate axes. In this case, the edge lengths of the typical cell of the 
rectangular Poisson line tessellation are independently exponentially distributed with mean 1; see [5]. 
Hence, we simulated 10 12 independent realizations of the random vector (X, Y) with X and Y i.i.d. 
Exponential(l). From the collection of cells obtained this way, we selected the 150 cells with the 
smallest area. These are shown in Figure [2] (left). Histograms of the deviation functionals a and r are 
shown in Figure [3j Both, the line segment shape of the cells with an accumulation around the origin 
and the peak at zero in the histograms for the deviation functionals are nicely visible. 

As discussed above, the area is only one measure of size of a cell. We investigate now the shape 
of small cells in the sense that their perimeter tends to zero. Therefore, the 150 cells with smallest 
perimeter were extracted from the sample generated above. These cells together with their deviation 
functionals are shown in Figure [2] (right) and Figure [4j respectively. In this situation, the maximum 
of the edge lengths also tends to zero, but Theorem [2] implies that a follows a uniform distribution on 
the interval [0, 1], which is also nicely visible in the histograms. The minimal and maximal sizes of 
area and perimeter of the cells included in the statistics described above are given in Table [TJ 

The Poisson line tessellations considered in this paper have natural analogues in higher space 
dimensions, the so-called Poisson cuboid tessellations for which we refer to [S]. Also for these tessel- 
lations the question about the shape of small cells can be asked. Natural candidates are the typical 
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Figure 3: Histograms for the deviation functionals of the 150 cells with smallest area. 





minimum 


maximum 


area 2D 
perimeter 2D 


1.79e-14 
1.06e-6 


8.46e-12 
3.52e-5 



Table 1: Minimum and maximum of the size of the cells considered in the statistics (planar case). 



cell of small volume, small surface area or small total edge length. Unfortunately, we were not able to 
derive the higher-dimensional pendants to Theorem [T] or Theorem [2] in full generality. This is mainly 
due to technical complications that arise for space dimensions d > 3 and cause that the Abelian- 
type theorem used in the proof of Theorem [T] can no more be applied. For this reason we carried 
out a simulation study concerning small cells in M 3 . For this purpose, 10 12 independent realizations 
of the random vector {X\,X2,X^) with Xx^X^^X^, i.i.d. Exponential(l) were generated representing 
the random edge lengths of the typical cell. From the sample of cells the 150 cells with smallest 
volume, surface area and total edge length were extracted. Histograms for the deviation functionals 
a = 3 min(Xi , X2 , X$) / {X\ + X2 + -X3) (again the factor 3 yields a value of 1 for a cube and the range 
[0, 1]) and t = max(Xi, X2, X3) for these cells are shown in Figure [5] The minimal and maximal sizes 
of the cells included in the statistics above are summarized in Table [2j 





minimum 


maximum 


volume 3D 
surface area 3D 
total edge length 3D 


3.20 e-15 
1.74e-8 
6.80e-4 


4.97e-13 
3.58e-7 
3.88e-3 



Table 2: Minimum and maximum of the size of the cells considered in the statistics (spatial case). 



5 



I 1 1 1 1 1 I 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 4: Histograms for the deviation functionals of the 150 cells with smallest perimeter. 



3 Proof of Theorem [T] 

3.1 Reduction 

We claim that without loss of generality we can restrict the proof of Theorem [T] to the case 7 = 2, 
q = 1/2 and Q = \{5l x + $L V ), where L x and L y are the two orthogonal coordinate axes. To show 
this, let us denote such a tessellation by PLT* and a tessellation with general parameters 7, q and Q 
by PLT(7, q, Q). We notice now that due to our assumptions on 7 and q there exists a non-degenerate 
linear transformation / = f{^f, q, Q) : M 2 — > M 2 such that / applied to PLT(7, q, Q) has the same 
distribution as PLT*. We notice further that the images under / of a parallelogram, a line segment 
and a point are again a parallelogram, a line segment and a point, respectively. Thus, the statement 
of Theorem [T] is invariant under non-degenerate linear transformations of the underlying Poisson line 
tessellation. This implies that it is sufficient to establish the statement for one special choice of 7, q 
and Q, viz. PLT*. 

3.2 Proof for PLT* 

Proof of ([5]). To simplify the calculations we work from now on with a = a/2 and translate the 
result afterwards to the original deviation functional a. To start with the calculation, we write the 
conditional probability as 

F(a > e\A < a) = y — - r — -. 
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Figure 5: Histograms for the deviation functionals a (first row) and r (second row) of the 150 cells with 
smallest volume (left), surface area (middle) and total edge length (right) in K 3 . 



In what follows, we consider the numerator and the denominator separately. The numerator can be 
written as 




dy&x 



The last expression can be determined by a straight forward integration procedure, which yields that 
asymptotically, as a — > 0, it behaves like a times a constant depending on e (the two logarithmic terms 



cancel out). For the denominator it follows that 

JJ e - x - y dydx = J J e- va(x+y> adydx 

A<a s 2 -t 2 <4 




oo 




/* 


—</as 


2 


oo 


a 






/• 


~ 2 


2 



(9) = l-(i±M^-|/^.(^)d, 

2 

Even if the integral in ([9]) looks rather innocent, its asymptotic behavior as a — > turns out to be not 
accessible with elementary methods as above. To overcome this difficulty we make use of a theorem 
of Abelian type. So, let F(s) := V s 2 + 4s, write 

oo oo 

1(a) -=^J e-^iVs 2 -4)ds= ^ e -^ s F(s) ds 

2 

and use now the Abelian Theorem [U Theorem 1.213] to conclude that 

(10) lim X( " } = 1. 



Putting together Q and ( 10 ) we see with a Taylor expansion of the exponential function that P(^4 < a) 
behaves asymptotically like \fa in the limit as a — > 0. Combining this with the asymptotic behavior 
of the numerator implies ([5]) for a and also for the original deviation functional a. □ 

Remark 3. At the beginning of the proof we used the very rough estimate e~ x ~ y < 1 for the integrand 
in the numerator. This can be refined considerably by explicit calculations. This is possible except of 
one integral whose asymptotic behavior can unfortunately not be discussed with an Abelian theorem. 
Moreover, all estimates we found did not improve the rate 0(a) as a — > 0. Progress in this direction 
would be necessary to optimize 

Proof of 0. At first, we re-write the conditional probability as 

5 (max{A, Y} > e, A < a) 



(r > e\A < a) 



P(A < a) 



The denominator is the same as in the proof of ^ and behaves like \fa as a — > 0. Next, we consider 
the numerator and write 

P(max{ X, Y} > e, A < a) < P(X > e,A<a) + ¥(Y > e,A<a) 

= 2P(A > e,A<a) = 2 J J eT x ~ y dydx. 

Since numerator and denominator both tend to as a — > 0, we may apply l'Hospital's rule. For the 
numerator we find 



d_ 
da 



poo roc ^—^—x 

/ / e~ x ~ v dydx = I dx -)• T(0, e) + lne < oo as a ->■ 0, 

Je JO Je x 
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where T( • , • ) stands for the lower incomplete T-function. Next, we turn to the denominator. Differ- 
entiation and rearrangement of terms yields 

7 rco r— -i poo 

-f- / I e- x - y dydx = -e- 2 ^ + -(l-Va)e~ 2 ^ e'^ s V ' s 2 + 4s ds 
da J Q J 2 J 



4 Jo 



--■T 2 



The term T\ has already been discussed in the proof of ^ above. Thus, we may concentrate on T 2 
and write 72(a) to indicate the dependence on a. Defining F(s) := sV s 2 + 4s we have that 



T 2 (a) = ^e~ 2 ^ / e-^ s F(s)ds 



4 



CO 







and can apply [H Theorem 1.213], which implies that ^2(0) behaves asymptotically as e~ 2 ^/a as 
a — > 0. Putting things together we see that the numerator, after differentiation, is bounded, whereas 
the denominator, again after differentiation, tends to —00 as a — > 0, implying that 

limP(r > e\A < a) = 0. 

a— >0 

This completes the proof of ([7]). □ 

Remark 4. Although the application of the Abelian theorem gives a precise asymptotic for the 
resulting integral, we are not able to determine from this a rate for P(r > e\A < a) in the sense of 
(§. This is because we have used l'Hospital's rule in the beginning of the proof and have considered 
the derivatives thereafter instead of the original expression as in the proof of ^ . 

4 Proof of Theorem [2] 

Applying the reduction step as in the proof of Theorem [T] we can again restrict our attention to 
the case of PLT*. Recall that for PIT* the random edge lengths X and Y of the typical cell are 
independent and identically distributed according to an exponential distribution with mean 1 so that 
the half perimeter length P = X + Y is Erlang distributed with parameters 2 and 1. In fact this was 
the reason for the factor 1/2 in the definition of P. Thus, 

(11) P(P < p) = 1 - (1 + p)e~ p for any p > 0. 

Without loss of generality assume that X > Y and have a closer look at the event {a > e, P < p}. 
If X € [0,p/2], then Y may range between 2^X and X, and if X E \p/2,p(l — s/2)] then Y ranges 
between ^^.X and p — X (the remaining case X > p(l — e/2) contradicts X + Y < p and Y > t^X). 
Thus, 

p/2 x p(l-e/2) p-x 

P(<r > e, P <p) = 2 J J e- x - y dydx + 2 J J e~ x - y dydx. 

^-x p/2 ^-x 



Evaluation of these integrals yields P(cr > e, P < p) = (1 — e)(l — (1 +p)e p ) and in view of (11) the 
exact distributional result 

P(cj > e\P <p) = 1 -e. 

Next, we notice that whenever p < e we have the following situation (without loss of generality we 
may assume that max(X, Y) = X): we want that X > e and X + Y<p<e&t the same time, 
which is impossible. Consequently, P(r > e, P < p) = holds for such p. Finally, this yields ^ and 
completes the proof. □ 



9 



Remark 5. We would like to point out that the first result of Theorem [2] has a well-known background. 
Namely, let X and Y be two independent and exponentially distributed random variables. Then X 
(or Y), given that X + Y = s for some fixed s > 0, is uniformly distributed on [0, s]. 
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